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1 Introduction 



The Bethe lattice, or the infinite Cayley tree, presents a hierarchical structure that greatlty simplifies 
some problems of statistical physics. It has therefore been widely used to obtain analytical results for 
problems that are otherwise intractable on Euclidean lattices. The physical relevance of these results 
is that the Bethe lattice is supposed to represent some mean field limit of Euclidean lattices of very 
large dimensions. 

The tree structure is also directly relevant to characterize the phase space hierarchical structure 
of some disordered systems like spin glasses jjj, where this structure plays an essential role in the 
understanding of the slow dynamics properties of these systems [g] |J Q. 

Random walks on Bethe lattices have already been studied either for their own interest || || , 
or in the context of polymer physics ||, where counting unentangled loops of a polymer in an array 
of obstacles is equivalent to counting random walks on Cayley trees. In this paper, we give the exact 
solution to this problem of random walk on the Bethe lattice and discuss the correspondance with 
Brownian motion on a space of constant negative curvature. 

The paper is organized as follows. In section ^, we consider a mapping between random walk 
on the Bethe lattice and a biased random walk on the half-line with reflection at the origin. We 
derive asymptotic expressions at large time through a simple argument, and write the exact explicit 
solution using old results of Mark Kac |1C]. Section |3] is devoted to biased Brownian motion on the 
half-line with a reflective barrier at the origin, which represents the continuous limit of the previous 
discrete model. Finally in Section [|, we discuss the relation between the Bethe lattice and hyperbolic 
geometry. We compare the properties of Brownian motion on a space of constant negative curvature 
to the results obtained in the previous sections. 



2 Random walk on the Bethe lattice and biased random walk on 
the half- chain 

Hughes and al. || have pointed out that a random walk on the Bethe lattice can be mapped onto a 
biased one-dimensional lattice walk with a reflecting barrier at the origin, provided one only considers 
the distance to the origin on the Bethe lattice as a function of time. In this section, we first rederive this 
correspondence to fix our notations. We then use the generating function formalism and a Tauberian 
theorem to obtain the large time asymptotic behavior of the probability of being at any given distance 
from the origin at time t. This result generalizes a previous one concerning the asymptotic behavior 
of the probability of returning to the origin at time t |] || 0. We finally use an old result of Kac 
[ |io| l to write the full solution of the problem of random walk on the Bethe lattice in an exact closed 
form. 



2.1 Master Equation 

We consider a random walk on the Bethe lattice of coordination number z > 2. We take the origin of 
the random walk as the origin of coordinates, and call generation n the set of all the z(z — l) n ~ 1 sites 
distant from the origin by n > 1 bonds (see Fig.l). 

The random walk is defined as follows. At each time step, the particle jumps with probability f-M 
to any of its z nearest neighbors. Therefore, the particle goes further from the origin with probability 
(^T^)j and goes closer to the origin with probability (j)- The corresponding Master Equation for 
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the probability f T {n) of being on any site of generation n after t time steps reads 

f T+1 (n) = (l - fr(n - 1) + ~/r(n + 1) for n > 2 (1) 

This has to be supplemented by specific rules at the boundary sites n = 1 and n = 

/ T+1 (l) = / T (0) + ~/ T (2) (2) 
/r+l(0) =~/r(l) (3) 

and by the intitial condition at r = 

f T =o{n) = S n>0 (4) 

oo 

The normalization reads /r(^) = 1 for any r. 

n=0 

The random walk on the Bethe lattice can therefore be described as an asymmetric one-dimensional 
random walk on the half-line n > with a reflecting barrier at the origin n = || . One may then use 
standard techniques for one-dimensional random walks. 

2.2 Asymptotic behavior at large time 

It is convenient to introduce the generating function 

oo 

F n (X) = y £X T f T (n) (5) 

T = 

The master equation (|l|) for f T (n) is then converted into the recurrence relation 

F n+1 (X) = jF n (X)-(z-l)F n ^(X) for n>2 (6) 

and the boundary rules (|2||3]) together with the initial condition (||) give the relations 

*2(A) = jFx(X) - zF (X) (7) 

F 1 (X) =£[f (A)-l] (8) 

Two linearly independant solutions of (|6|) are [r(A)] n and [r'(A)]' 1 where r(A) and r'(X) are the solutions 
of the characteristic equation 

r 2 - + (z - 1) = (9) 
A 

As we must only keep the regular solution in the limit A — > (||), we get for n > 2 



^(A) = [r(A)] n - 1 F 1 (A) with r(A) = — 



l-\/l-4(^-lM 2 



(10) 



The two first terms -Fi(A) and -FolA) are then determined by (f 
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Fo{X) =-, r y ' (11) 

1_ /l_ 4( 2 _1UJ \ 

^i(A) =t \ = (12) 



The probabilities f T (n) may then in principle be obtained by expanding the previous expressions (10- 
^) in powers of A (||). This has been done for the function -Fo(A) to obtain / T (0) in an explicit form 
involving hyper geometric functions || |7]]. This has not been done for the general case n > 0, where 
it is preferable to use some other method to get an exact explicit expression of f T (n) for any n (See 
next paragraph |2.3| ) . 

However, the asymptotic behavior at large time of f T {n) for any fixed n may be easily found from 

the generating function using a Tauberian theorem. 

z - 2 2 

Introducing the notations v = and a = —yz — 1, the generating function reads for n > 1 

z z 



F n (X) = 2 
with the auxiliary function 

9n(x) = 2 

that admits the power series expansion 

9n{x 



1 



v + \J\ — (c*Ap 



za 
2x~ 



(i-vT 



{fc>0;(fe-n)even} 



(13) 



(14) 



(15) 



Since we have the simple relation f T (n) = c T (n)a T , the asymptotic behavior of the probability 
f T {n) at large time r is directly related to the behavior of the coefficients c^(n) for large order k. The 
latter can readily be obtained from the expansion of g n (E3) near x = 1 



9n{e- E ) 



2 (za\ n 



1 - ( n + - ) V2E + 0(E) 



Using the power series expression (|i5|), we get 

X; c fc (n)(l-e- fcS 
{fc>0;(fc-n)even} 



1\ f za 



~ n + - Hr V2£ + 0(£) 



(16) 



(17) 



The non-analyticity in \/E implies that the coefficients c^(n) present the following algebraic decay at 
large order k for (k — n) even 



Cfc(n) ~ 



A(n) 



fc ^°° ki 



The prefactor yl(n) can also be obtained from ( |17|) by converting the sum into an integral 



1 \ / za V 



A(n) = - (n + - — 



3 
22 



2 foo 



du 



l-e~ u 



22 (1 + nu) fzaV 

7^ v 2 vtJ 



We finally get for (k — n) even 



(18) 
(19) 

(20) 



a' 



f T (n) = c T (n) a T ~ A(n) — 



(21) 



4 



or more explicitly 



2§ 

fr(n) ~ -p 



T— >oo ./7T V 2 — 2 



1 + n- 



z-2 



s/ Z — l) — r e 



-t In 



3 
T2 



(22) 



This formula valid for n > 1 generalizes the result obtained previously for the case n = 



2§ z(z-l) 1 " rln 
/ T (0) ~ -= / ^ — e 



r— x v 77 iZ - 2) 2 r f 



2VJ^1 



(23) 



2.3 Exact expression for f T (ri) 

In [10], Kac studied the asymmetric random walk ([l]) with the reflective boundary condition (§-§), 
but in the domain of negative drift 1 < z < 2. Using our notations, the explicit expression obtained 
in this paper for the initial condition (Eh reads 



{0<z<2) 



(») 



2 

+- 

7T 



(2-z) 
z-1 



5n,0 + z(z-l) n 1 (1 - 5 n>0 ) 



l+(-l)" +T 
2 



5 nfl + (z - 1) 2 (1 - <$„ „) 



'2v^T\ r 



I(n,r) 



where 

/(„,r) = [l + (-l) n+ H 



d0 (cos 0) T 



tarn 



(£_2) + tan 2 # 



cos no + 



z — 2\ sinn6 



Sill! 



cos 9 



(24) 
(25) 

(26) 



The first term (|24|) can be easily recovered as the normalized "stationary" (invariant under the time 
translation r — > r + 2) solution of (|l|) in the case of negative drift 1 < z < 2. 

For the case of positive drift z > 2, which we consider in this paper, there is no normalizable 
stationary solution. The exact expression for f T describing the random walk on the Bethe lattice 
therefore only contains the second term (p5|) 



(*>2). 



n) 



z-1 



S n ,o + (z - 1) 2 (1 - «5„,o) 



I(n,r) 



(27) 



For the particular case n = 0, the integral 1(0, r) fl26|) may be computed through the change of 
variable x = cos 2 # 



1(0, 2t) = 2 / d# (cos 0) 



s2t 



tan 2 6> 



(*=2) + tan 2 fl 



T-l 



1 — x)i 



which is a standard integral representation of the hypergeometric function 



1 3 



1(0, 2t)=B[t + -,-) F 1, r + -, r + 2 



2' 2 



1 



4(z - 1) 



We therefore recover the expression given in ||| Q 



(0 ) 



1\ / v 7 ^^' 



2r 



r(2r + l) 



r(r+l) T(r + 2) 



1 4(2-1) 
F l, r +-,r + 2, 1 9 ; 



(28) 



(29) 



(30) 
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We may also recover from the expression fl27|) the asymptotic expression at large time given previ- 
ously in (p2|). Indeed for large r, the integral I(n, r) (26) is dominated by the vicinity of 9 = 0. The 
estimation of the leading order 



I(n,r) 



[1 + (-l) n+r ] / d6 e~ T - 9 2 
Jo 







1 + 


(^H 







l+(-l 



,n+T 



z-2 



1 + 



z-2 
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3 

T2 



(31) 
(32) 



gives back (p2[). 

3 Biased Brownian motion on the half-line 

We consider the continuum counterpart of the biased random walk described in section (|l|). The 
probability density P(x, t) for a Brownian particle of diffusion constant D submitted to a constant 
positive drift \\i satisfies the Fokker-Planck equation 



dt dx V dx 
It must be supplemented by the reflection condition at x = 

dP 



D—-fiP)(x = 0,t) = 



(33) 



(34) 



to ensure conservation of probability, and by the initial condition 

P(x,t ^0+) = S(x) (35) 

3.1 Expression of the probability density P(x,t) 
It is convenient to perform the transformation 

P(x,t) = ip(x,t) (36) 

in order to cast the Fokker-Planck equation ( |33"1) into the Euclidean Schrodinger equation 

dip 



dt 



-Hip 



The corresponding Hamiltonian is a free one up to a constant shift 

j2 ,,2 



H 



dx 2 AD 



The boundary condition and the initial condition read respectively 

(x = 0,t) = 



dip 

~dx ~ 2~D 



and 



ip(x,t^ 0+) = 5(x) 



(37) 

(38) 

(39) 
(40) 
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Using plane waves, it is easy to construct an orthonormal basis of eigenvectors {^(x), k £ [0, +00 [} 

,2 \ 



Hip k = ( Dk 2 + I ip k 



(41) 



satisfying the boundary condition 



1 



—ikx 



/i + 2iL>/c eikx 



'2ir V - 2iDfc 

The Green function tp(x,t) solution of ([37-39-4C) may be expanded onto this orthonormal basis 



(42) 



ip(x, t) =< x\e 



-tH 



>= 



dk tp k (x) V'fe(O) e 



(43) 



As a side remark, let us mention that in the case of negative drift < 0, in addition to the 
continuous spectrum (f42|), there is also a bound state solution of zero energy which must be added to 
the expansion ([43]), This zero-energy state corresponds to the aforementioned stationary solution of 
Kac the discret case (|24|). 

However, in the case of positive drift /j, > that we consider here, there is no stationary normal- 
izable solution, and the spectrum in purely continuous (f43|). The resulting integral 



1 



ip(x, t) - — I dk 
can be computed through the following trick 



' h 2 — ik-tL\ 



2 

417 2 



(44) 



ip{x,t) 



e 4D " I(x) 

{ dx 2 2D8x { ' 



(45) 



The integral I(x), being the Fourier transform of the product of a Gaussian by a Lorentzian, may be 
written ctS Si convolution 



I(x) = - dk e ikx 
vr J -00 



1 



k 2 + J±L 



D 



i 

dy e 2D 



\x - y\ 



1 



Going back to the probability density (|36|), we finally obtain 

,2 



4L>t 



(46) 



P(x,t) 



V^Dt 



x- 



e IDt-JLl dye 2D [V X) e 4Dt 
2D J x 



(47) 



It is useful to transform the previous expression into 

,2 



P(x,t) 



— t a 

e 4D —x 
* e2D 



2y/^(Dt)l 



-My 

dy y e 2D 



ADt 



(48) 
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The asymptotic behavior at large time t immediately follows 



P(x,t) ~ — r eD / dyy e 2ZT (49) 



2 



7T /x 2 V 2i?y t i 

Comparison with the corresponding result for the discrete case (|2 



2^{Dt) 

2 



/ T (n) * — 1= — — — ) [1+n— e n ^-^e (51) 



72 

shows that even if the two expressions behave qualitatively the same as t~2 e _j4 * as a function of time, 
and as (1 + nC) e nB as a function of space, there is no direct detailed correspondance at this stage. 
If one naively tries to identify x = n and t = r, then there is no way to find two functions fj,(z) and 
D(z) to make equations ( |50|) and (]5l| ) identical. The relation between the random walk on the Bethe 
lattice and the continuous model (|33|) studied in this section therefore needs further clarification. 



3.2 Discussion of the continuous limit 

To define properly the continuous limit of the random walk on the Bethe lattice, we must introduce 
some unit length Ax and some unit time At into the Master equation (|l]) 

f(x, t + At) = (l- f(x - Ax, t) + ~/(x + Ax, t) (52) 

The Taylor expansion 

dj_ = ( {Axf \ cP£ _ f z-2 Ax\ df_ 

dt \ 2At J dx 2 V z At) dx [ 1 

gives the Fokker-Planck equation (|3^) in the limit 

Ax -> 0+ ; At -» 0+ ; z -► 2+ (54) 

with 

{ -^ = D and Z -^^ = , (55) 

2At 2 At ^ y J 

where D and \i are fixed. This continuous limit therefore requires an analytic continuation to non- 
integer coordination number z in order to take the limit z — > 2 + . With the prescription ( |54"f]5"5| ) , the 
result (|22|) for the discrete case indeed corresponds to the result ( |50t ) of the continuous case. 

Let us now write for the discrete case the analogue of the transformations that we performed for 
the continuous case in order to establish some correspondance with the work of Clark et al. [11]. The 
analogue of the transformation (36) to make the walk (|l|) symmetric 

/r(n) = e nln v / ^T f r(n) (56) 
and the extraction of the shift factor analogue to the shift present in (j3§|) 

-r In 



f T (n) = e W*-l/ p T (n) (57) 
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transform the drifted random walk with reflection at the origin into a symmetric random walk 

with a partial absorption at the origin 

1 1 

p T+ i(n) = - p T (n - 1) + - p T (n + 1) for n > 2 (58) 
p T+1 (l) =ipM + \p T {2) (59) 

Pr+l(0) =\prO) (60) 

where 7 = 2(2-1) denotes the reflection coefficient at the origin. The continuous limit of this random 
walk with partial absorption at the origin is indeed a quantum mechanical free problem on the half-line 



with some mixed boundary condition (|39|), The detailed discussion of this point contained in [If] is 



equivalent to our approach (see |5J-| 



The meaning of the continuous model ( p3[ ) to represent the random walk on the Bethe lattice is 
now clear. We will now present another continuous model where the effective drift comes from the 
geometry of the underlying space itself. 

4 Bethe lattices and Hyperbolic geometry 



A link between Cayley trees and hyperbolic geometry has already been emphasized in [12] |13|] through 
the introduction of tessalations. Indeed, for any integer z > 3, one can construct a tessalation of the 
Poincare upper half-plane with polygons with z sides. Joining the centers of adjacent polygons gives 
a Bethe lattice of coordination number z. In the following, we will not use this approach, but rather 
establish some intrinsic link between Bethe lattices and spaces of constant negative curvature without 
introducing any tessalation. 

4.1 Hyperbolic geometry 

An N-dimensional Riemannian manifold of constant negative Gaussian curvature K = — may be 
described by the metric 

ds 2 N = dr 2 + a 2 sinh 2 - da 2 N _ 1 (61) 
a 

where r G [0, +00] measures the distance to the origin, and where da 2 N _ 1 denotes the metric of the 
unit-sphere Sw-i. For example, da 2 = d6 2 is the metric of the unit circle in terms of polar angle 9, 
and da 2 = d6 2 + sin 9 2 d(j) 2 is the metric of the unit sphere S2 in terms of the spherical angles (#,</>)• 
The volume element dV is covariantly defined as 

dVjy = (a sinh— J dr (iiljv-i (62) 

where d£lN-i is the surface element of the unit-sphere Sm-i\ for example d£l\ = dO and d£l2 = 
sin OdOdcf). 

In particular, the volume of a ball of radius R reads 

V N (R) = ftjv-i a N I* dx (sinhx) 7 ^ 1 (63) 
Jo 

where £In-i is the total surface of the unit-sphere Sn—i', for example S7i = 2tt and ^2 = 47r. This 
volume has the important property to grow exponentially with the distance R 

R-' N - V 



V S (R) « e \ <* J (64) 

KS>a 
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in contrast with the power dependance R N in Euclidean space of dimension N. 
In the same way, the number of sites on the Bethe lattice up to generation n 



N(n) = l + z + z(z-l)+z{z-lf + --- + z{z- I)™" 1 = { J (65) 

grows exponentially with the number n of generations 

M(n) oc e ^M^-l) ( 66 ) 

This is why the Bethe lattice is often considered to be like a Euclidean lattice of infinite dimension. But 
it is certainly more interesting to relate it to hyperbolic geometry which presents the same property 
(|64|). More precisely, we may introduce some lattice spacing Ar on the Cayley tree so that the sites 
of generation n are at distance r = nAr from the origin. Then the exponential dependence of d66| ) 
corresponds to ([64] ) in the following continuous limit of the Bethe lattice 

Ar -► + ; n -> oo ; z -> 2 + (67) 

with 

N-l , \ 

nAr = r and = (68) 

Ar a v ' 

where iV and a are fixed. 



4.2 Hyperbolic Brownian motion 



The radial part A r of the Laplace operator on the N-dimensional Riemannian manifold of constant 
negative Gaussian curvature defined by the metric (|6l|) reads 







(sinh^- 1 Or 



sinh • 



N-l 



0_ 

Or 



(69) 



On this manifold, free Brownian motion starting from the origin is defined by the diffusion equation 
for the Green's function Gjv(r, t) 

&Gn DA r G N 



and the initial condition 



Of 



G N (r,t) 5(r)- 



1 



t->o+ r Ojv_i 

The normalization of the Green function G^(r,t) then reads for any time t 

N-l 



/r+oo / 
dS G#(r, t) = Qn-i J dr fa sinh 

The solution reads respectively in two and three dimensions 

_ Dt 

G 2 (r,t) 



G N (r, t) 



y 



4a 2 



4V2a(iTDt) 



3 

2 Jr 



°° ye ^Dt 

dy 



cosh (|) - cosh (I) 



(70) 
(71) 

(72) 



(73) 
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and 



G 3 (r,t) 



Dt 



r e 



4Dt 



3 7" 

la(irDt)2 sinh- 
a 



(74) 



Consider now the probability density Pjv(r, t) to be at time t at a distance r from the origin 

N-i 



PAr) = Sn-i a sinh 



normalized with respect to the flat measure dr 



G N (r,t) 



(75) 



dr P N (r,t) 



This probability density satisfies the Fokker-Planck equation (17 



dP _ d 
dt dr 



and the initial condition at t = (|7T| ) 



or 



P N (r,t) t -j + 5(r) 



coth ( - ) P 
a 



The solution reads respectively in two and three dimensions (|73|-74) 



2V2vr(-Dt) 



sinh 



y e 



4L>i 



cosh (|) - cosh (I) 



P 3 (r,t) = a- 



Dt 



r sinh — e 4Dt 



(76) 



(77) 



(78) 



(79) 



(80) 



2y/Tr(Dt)% a 

At large distance from the origin r 3> a, the Fokker-Planck equation (|77| ) corresponds to a one- 
dimensional diffusion with a constant drift \x = D N ~^ (H; 



<9P 9 
<9t <9r 



iV- 1 



P 



(81) 



One may check that this identification of the effective radial constant drift fi = D Z 1 at large 
distance on the hyperbolic space is entirely consistent with the two continuous limits of the Bethe 
lattice previously defined in ( |54"f]5"5| ) and (67-58). 

The solutions (79|-80|) in two and three dimensions read approximatively at large distance 



Dt 

P2(r,t) ~ r e2a 

r » a A^(Dt)l 



y-r 



V 



dy y e 2a e 4Dt 



Dt 

>r»a 4^(^)1 



a e ±Dt 



r e« e 



(82) 



(83) 
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to be compared with (^). These expressions at large time 



P 2 {r > a, t —* oo) ~ 



3 r e2a -g. e 



1 

is 



4^2 



(84) 



P 3 (r,t) ~ 

r '» a 



r Z3i 

- 1 2 

T r ea e a 



t2 



are the exact analog of expression (|50|) at large distance 



(85) 



2D 

P(x > ,t -> oo) = 

fi 2/iVvrLI 



*5 



(86) 



for /i = D N a l . This confirms the equivalence with biased Brownian motion on the half-line. 



5 Conclusion 

The exact solution (|27| ) to the homogenous random walk on the Bethe lattice was obtained through 
a mapping onto a biased one-dimensional random walk on the half-line that describes the "radial" 
dynamics on the tree. This mapping to a one-dimensional lattice can of course be extended to inho- 
mogenous random walks on the Bethe lattice that still preserve the radial invariance on the tree |j] 
[ |j~4| l |nf |, but does not hold any longer as soon as two sites belonging to the same generation are no 
more equivalent. This is in particular the case when one considers disordered problems on the Bethe 
lattice. This is why the relation with hyperbolic geometry that we presented in section ^| is much 
more profound. On a Riemanian manifold of constant negative curvature, we have at our disposal 
not only the radial coordinate that corresponds to the generation number on the Bethe lattice, but 
also angular coordinates, that correspond to the "degree of freedom" inside a given generation of the 
tree. In this paper, we have discussed in detail the radial correspondance, but it would certainly be 
interesting to consider also more precisely the angular part, and to study in particular the continuous 
limit of some disordered systems on the Bethe lattice. 
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